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Abstract 

We provide a compendium of techniques that can be used to compute 
0{Nfa'^) corrections to low-energy electroweak processes. Specifically, these 
are the 2-loop electroweak corrections containing a light fermion loop. It is 
shown that the vast majority of such corrections can be reduced to expressions 
involving a universal master integral for which an exact analytic form in dimen- 
sional regularization is given. The only exceptions are certain photon vacuum 
polarization diagrams. Examples are presented for diagrams that occur in a 
variety of processes of practical interest. 



1 Introduction 



Among the most precisely measured electroweak observables the electromagnetic cou- 
pling constant, a, the muon decay constant, G^, the anomalous magnetic moment 
g — 2, and the ratio of charged to neutral current cross-sections in semi-leptonic neu- 
trino scattering, feature prominently. Because of their high experimental precision 
a and are commonly used as input in calculations of electroweak radiative cor- 
rections. Until relatively recently the experimental errors on a and were entirely 
negligible compared to those on the third input quantity, Mz- Now, however, both 
Gfj^ and Mz are both measured to 2 parts in 10^. Any further improvement on the 
accuracy of Mz would mean that the measurement of the muon lifetime in principle 
limits the accuracy with which electroweak theory can be tested. 

In this arena, theoretical calculations currently lag experiment. Calculations are 
not available that could exploit the 2 x 10~^ accuracy available on and Mz- 
It is therefore of technical interest and practical necessity to develop electroweak 
calculations that approach the 2 x 10~^ accuracy of experiment. While a full 2-loop 
calculation would be ideal, it is likely that the analytic form of such results would 
be too complicated to usefully write down. One would therefore probably wish to 
rely on computer algebra and integration to obtain numerical results [0]. However 
it is still true that certain dominant subclasses of 2-loop diagrams can be calculated 
and written down in compact analytic form Such results are, of course, always 

preferable to numerical ones. 

The set of 2-loop Feynman diagrams containing one fermion loop are such a class. 
These we will refer to as (9(A^/a;^) corrections where Nf is the number of fermions 
and Ng = Nf/8 is the number of complete generations. In these corrections we will 
assume that fermions are massless which is an excellent approximation for all but the 
third generation. The 0{Nfa'^) results can be used either with Ng = 3 and mass 
corrections computed for the third generation or Ng = 2 and corrections for the third 
generation computed separately. 

Because Nj is quite large — up to 24 — the 0{Nfa^) corrections are expected 
to be a dominant subset of 2-loop graphs and since Nf provides a unique tag, the 
complete set of Oi^NfO^) corrections contributing to a particular physical process will 
form a gauge- invariant set. Diagrams of this type have been considered in connection 
with the muon anomalous magnetic moment 0. 

A priori the 0{Nfa^) corrections can be expected to contribute at the level of 1.5 x 
10~^. Without their calculation and inclusion, theoretical predictions are uncertain 
at this level. 

To test the Standard Model, one needs at least four well-measured electroweak 
observables. Three of these are used as input to the model and the fourth is then 
predicted and can be compared with experimental determinations. In addition, re- 
lating electroweak parameters, such as Mz-, extracted at high energies with those, 
such as a and G^, extracted at low energies inevitably involves an uncertainty aris- 
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ing from the hadronic contribution to the photon vacuum polarization, ll'^^{0). This 
hadronic uncertainty can be reduced by improved measurements of the cross-section 
cr(e+e~ hadrons) and the use of dispersion relations or eliminated completely 
by sacrificing one high precision observable [|10|] . 

At present the fourth well-measured electroweak observable, needed to test the 
Standard Model at very high precision, is missing but it is not unlikely that quantities, 
such as the polarization asymmetry. Am, or — 2 of the muon, could fulfill this role 
in the foreseeable future. In any event the calculation of the 0(iVja^) electroweak 
radiative corrections are a necessary prerequisite to be able exploit the extremely 
accurate experimental measurement that are the legacy of LEP. 

2 The Master Integral 

In what follows the general coupling of a fermion to a vector boson will be denoted 



where jl and are the usual left- and right-handed helicity projection operators 
and Pl and are the corresponding coupling constants. We assume throughout an 
anticommuting 75. A Euclidean metric with the square of time-like momenta being 
negative will be used and all calculations will be done in R^=i gauge. The sine and 
cosine of the weak mixing angle 9w will be denoted se and ce respectively. 

We will be interested in the 0{Nfa'^) corrections for massless fermions. In cor- 
responding Feynman diagrams the fermions couple only to vector bosons and never 
to Goldstones or the physical Higgs particle. Many, but by no means all, such dia- 
grams are obtained simply by inserting a fermion loop into the boson propagator of 
a one- loop diagram, see Fig. la. Because the original one-loop diagram is logarith- 
mically divergent, the fermion loop insertion will need to be calculated to 0{n — 4) 
in dimensional regularization, where n is the dimension of space-time. It turns out, 
however, that for all OlNja"^) diagrams for low-energy processes, it is possible to 
obtain expressions that are exact in n. 

For the massless fermion loop insertion, it may be shown that 





T 

T{n) 



J [n-l)J i'n''q\q + pY 



Y2l n 



= - {P%. - P.P.) (PA + (3nf3'nH{7rp')^-'^^T (2 - ^)(2) 
where Pl, Pr and of the couplings of the attached vector bosons. 
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A number of special cases of low energy 2-loop integrals have appeared in the lit- 
erature |[TT|-[13|. By using the projection operator techniques described in the follow- 



ing sections, most C(A^/a^) diagrams occurring in low-energy electroweak processes, 
including box diagrams, can be reduced to expressions involving a general master 
integral of the form 

lii klmnM')= I ^> ^ ! ^ f3) 

provided the diagram contains at least one massive boson. The only exceptions to 
this are certain contributions to the photon vacuum polarization. In that case some 
diagrams must be calculated retaining the fermion mass, m/, and then applying the 
asymptotic expansion of ref. |]13[ before taking the limit — > 0. Other 0{Nfa'^) 



diagrams that cannot be expressed in terms of the integral are of the pure QED 
type containing only internal photons and fermions. The relevant master integrals 



for this case can be found in ref. 12 



In eq.(^, the integration over q can be performed using standard Feynman pa- 
rameter techniques and yields 



d-q 1 TTi-^ r(/ + m-|jr(5-/jrf§-m 



^2 [g2]«[(g + p)2]m [p2]«+m-f T{l)T{m)T{n - I - m) ' 

The resulting integral with respect to p in eq.(|^) is independent of angle and hence 

1 _ 27r?-2 p2n-2j-2l-2m-l 

m2 [p2]W™-f [p2 + M2]fc ~ P^nJ X ^ [p2 + M2]fc 

vrf ~2 Y(^n — j — I — m)T(k + j + l + m — n) 

(5) 



from which it follows 
/(j, k, /, m, n, 



X 



(^]\^2'^k+j+l+m-n 

r{n-j -I- m)r{k + j + / + m - n)r + m - § j r (^^ - r - 

r r(A;)r(/)r(m)r(n -l-m) 

(6) 



3 Self-energy Diagrams 

For electroweak processes, such as muon decay, in which the external fermions can 
be considered massless, only the self-energy diagrams of vector bosons need to be 
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considered. General techniques are available to reduce the tensor integrals that occur 
to standard form factors for which scalar integral representations are known. For 



the present case of zero external momentum, q = 0, the vector boson self-energies, 
n^i/(Q'^), can only take the form 

n^.(o) = 5^,F 

where F is a function of the internal masses only and may be obtained from the tensor 
integral representation of ll^^{0) by means of the projection operator, Sfj^^/n. Thus 



n 



^A^)- (7) 



The resulting scalar integral can always written in terms the master integral, /(j, k, /, m, n, M^), 
of eq.(|3) and results obtained that are exact for all n. 



4 Vertex Diagrams 

Again for processes with massless external fermions the only relevant vertex correc- 
tions are those involving vector bosons and these will necessarily be purely vector and 
axial-vector in character. A general vertex correction can then be represented as 

where Vl and Vr are functions only of the internal masses. The tensor integral 
representation of can easily be obtained by standard techniques and from it the 
scalar integral representations of Vl and Vr follow by means of projection operators. 
Thus 

V^L,R = -^Tr{V;7^7^,i}. (8) 

where Tr{7^7^} = 4n is assumed. This method for directly obtaining the scalar 
integral representation of the vertex form factors is particularly convenient when 
computer algebra is being employed. Once again the resulting scalar integrals can be 
written in terms of the master integral, /(j, fc, /, m, n, M^). 

4.1 Example 

The diagram shown in Fig.|I](a) occurs in the calculation of the 0{Nfa'^) electro- 
magnetic charge renormalization. It will be assumed that the external fermion has 
^3 = +1/2 in which case the fermion to which the photon couples has = —1/2. 
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(c) 



Figure 1: Feynman diagrams treated in the examples. 



We will denote the charge of the external fermion by Q. After obtaining the tensor 
integral for the Feynman diagram by usual methods, applying the projection operator 
as in eq.(^) gives the scalar integral representation 

-2 \ 2 
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Qse- 



n-2) [ dJ'p [ dJ'q [3p^q^ - 2{p ■ qY + q\p ■ q)] 



2 N 2 



1^2 p2^p2 ^ M^Yqi^q^^y 
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Qse 



in-2) 
2n 



X < -1/(1, 2, 2, -1, n, M^) - ^/(-l, 2, 2, 1, n, M^) + ^J(0, 2, 1, 1, n, M^) 



1-11^119 



(l|^)1«-r(4-.)r(2-|)r(|-: 



exactly. These diagrams and their related topologies do not vanish when summed 
over a complete fermion generation. 



5 Box Diagrams 

For 4-fermion processes one-loop box diagrams are finite but at 0{Nfa'^) they de- 
velop logarithmic divergences. A useful set of identities for calculating one-loop box 
diagrams appears in ref. |T^. They are, however, valid only for n = 4 because of their 
intended use at one-loop. For general n it may be shown that these relations become 

[7p7M7a7L,R]i[7p7;.7<77L,i?]2 = 45^;^[7a7L,i?]i[7a7L,i?]2 + - 4)[7p7L,R]i[7i^7L,i?]2 (9) 

[7p7M7a7L,/?]l[7p7i.7<77i?,L]2 = ^bvlLAlbf^^RA2+ - ^)bf^^L,R]lbvlR,L]2 (^0) 
[7p7p7a7L,/?]l[7<T7z^7p7L,R]2 = ^bvlLAlbt^lLA2+ i.^ - ^)bt^lL,R]lbulLA2 (H) 
[7p7M7a7L,i?]l[7<77;.7p7R,L]2 = 45^;^[7a7L,i?]l[7a7iJ,L]2 + («-4)[7p7L,/?]l[7i^7i?,L]2(12) 

where the square brackets [ ]i and [ ]2 indicate that the enclosed 7- matrices are 
associated with the external fermion currents Ji and J2 respectively. 

The general 0{Nfa'^) box diagram for massless external fermions at g = there- 
fore takes the form 

Bi_wJlfiJ2u = B ■ JlaJ2a- (13) 
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Here -B^i, is a tensor integral and the product of JiaJia can be constructed from one 
or a combination of the 7-matrices appearing in eq.(|9|)-(|T^). As for the self-energy 
contributions, the tensor integral B^j^y, can only be proportional to 5fj,u and hence, 
using the projection operator method of section 2, the scalar integral representation 
for the form factor B is seen to be i? = {6^u/n)B^u = B^^/n. In this case the 
identities (p|)-(|T2|) simplify to become 



6 y (5n - 4) 

-^[1 pl 1,1 alL Alb plyl'rlLA2 = [7a7L,/?] 1 [T^TL,/?] 2 (14) 

^[7p7M7<T7L,/?]i[7p7i.7<T7i?,L]2 = [7a7L,i?]i[7a7iJ,L]2 (15) 



n 



[7p7M7<77L,H]i[7a7!.7p7L,H]2 = [7a7L,R]i[7a7L,ij]2 (16) 



^[7p7M7<77L,/?]i[7a7v7p7/?,L]2 = ^^^^ '^\ lalLAibclRA2- (17) 
5.1 Example 

The diagram shown in Fig.0(b) occurs in the 0{NfQ^) corrections to muon decay. 
Using (1) and (|T5|) the resulting the Feynman diagram can be written as 

f g'^ Y r . r {n - 2) f f d^'q 

bc.lL]JlalLi 



9' f 9' \\ 1 r 1 (^-2) 



(M| - M^y 



i^^|/(0, 1, 1, 1, n, M|) - /(0, 1, 1, 1, n, M, 
+ (M|-M^)/(0,2,l,l,n,M^ 



g' f g' 



[7a7L]J7.7L],(vrM^)-^T(e) 



x<i(s^ + lnc2)h+ef^-7))+eln'cf 



up to 0(1) in e = 2-n/2. 



6 Derivatives of Self-energy Diagrams 

The derivatives of self-energies with respect to the square of the external momentum, 
k"^, are required for the calculation physical processes. They are often thought of as 
wavefunction renormalization factors although they occur even when wavefunction 
renormalization has not been explicitly performed [|T^. As we are concerned with 
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low-energy processes and massless fermions, these derivatives will only be required 
at k"^ = 0. In that case the scalar integral representation may be obtained by the 
following method. 

For a vector boson the self-energy may be written in the form 

n,.(fc^) = - ^) Urie) + n,ik^) (18) 

with Ut and being the transverse and longitudinal form-factors respectively. 

^Tik^) = - ^1 nM^(^') (19) 



k^kij 



^L{k') = {^]n,^{k'). (20) 



For massless fermions self-energy diagrams take the form 

S(fc) = {FL{k^)^L + FR{e)^n) (21) 

and 

FLAk') = -^Tr{E{k)^^n,L} (22) 

In either case given expressions for the self-energies in the form of tensor integrals, 
the scalar integral representations for the corresponding form factors fl^,, IIt or Fl, 

are easily obtained. In general, then, these scalar integral representations may be 
written 

F(P] ^ f f ^ F{k'^,p'^,q'^,k-p,k- q,p- q) 

(23) 



Z7r2 J ivr2 [(A; + p)2 + M2]^i[(A; + g)2 + M|]^2[(fc + p + g)2 + M3 



where Ni, N2 and are integer powers. 

In the region of small k"^ we may expand the propagators as a power series in {k -p) 
and {k ■ q). Thus, for example. 



[{k + pf + M2]^2 [A;2 + p2 + JVf2]JV2 ^ ^^2+^2 + ^2 



[{k + p + qf + M|]^3 [A;2 + (p + q)2 + V fc^ + + g)2 + m. 



with the result that the integral in eq.(p3|) becomes 
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Upon integration terms in which i + j is odd vanish. Terms in the integrand with 
i + j = 2 are {k ■ vi){k ■ V2)F{k'^, q^.p^) in which vi and V2 are either p or q. Because 
F is a function only of squared momenta the integration of p^quFiji^, q^iP^) can only 
yield a result proportional to 5^,^. Its coefficient may be extracted using the projection 
operator 5^u/n as in section 2. It follows then that the terms {k-Vi){k-V2)F{k'^ , q^.V^) 
in the integrand may be replaced by k'^{vi ■ V2)/nF{k'^,q'^,p'^) without changing the 
integral. 

In a similar way for i + j = 4 one can make the replacement 
{k ■ vi){k ■ V2){k ■ vs){k ■ Vi)F{k\ q^p") 

,^4 (^1 • ^2)(f3 • ^^4) + (fl • VZ){V2 ■ Va) + {Vl ■ Va){v2 " Vz) ^f.^j 2 2n 
^ 7 ^ y1 iP ) 

n[n + 2) 

where again the Vi are p or q. 

The integral is thus transformed into a form whose only k dependence is enters 
only through k^. 

Differentiation of the integrand with respect to fc^ can now be performed and 
setting k^ = yields a result expressible in terms of the master integral (^. 

For most practical purposes z + j = 4 is as high as one needs to go in the series 
expansion. The general term in this series can be found in fl^ and requires repeated 
application of the momentum-space d'Alembertian operator. 



6.1 Example 

A self-energy diagram for a massless fermion is shown in Fig.lc and may be written 

After obtaining the scalar integral representation for A^lk"^) and expanding the prop- 
agator as in ( p^ up to {k ■ pY one obtains 

1 / ^2 ^n-2) [ dy f (T^q 1 



Ak:^ \IQtx^) {n-l)J in^ J in'^ [p'^ + M^^Yq^q + pf 

~^ „5/L2 , y'^ ' P) 



AL{k' 



/^2 _|_ p2 p2(^/^2 _|_ p2^3 

After the replacement {k -pY — > /c^p^/n the derivative with respect to /c^ at /c^ = is 
found to be 

dAL{k^) 



/c2=0 



c/2 Y {n-2){n-A) f (T^p f d'^q 1 



167r2/ 4n J m"^ J m"^ ^(p^ + M'^Y q^ {q + p)"" 

V (n-2)(n-A) , , , 

(l|^)^'<-r-^r(5-„)r(2-^)rQ-i 
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exactly. This result is related to vertex diagrams by Ward identities that serve as a 
useful check. 
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